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_i , Abstract. We present form factors for a wide list of integrable models which in- 

cludes marginal perturbations of SU (2) WZNZ model for arbitrary central charge 
and the principal chiral field model. The interesting structure of these form factors 
is discussed. 
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Q\ ' 1. Introduction 

m ; 

The form factor bootstrap is powerful method of study of integrable field theories. 
There are many models for which the form factors are known, the most important 
examples are Sine-Gordon model, its asymptotically free limit: SU(2)- invariant 
Thirring model, 0(3) - nonlinear sigma model, ST/ (A^) -invariant Thirring model 
[1]. This list seems to be representative enough. However, for better understanding 
of integrable field theory we still need some more examples. Asymptotically free 
theories arc of primary importance for further investigations in the field. They 
allow important type of quantum symmetry: Yangian symmetry [2], which is in less 
modern language the same as Lusher nonlocal charges with interesting classical limit 
[3] . The knowledge of exact out of shell solutions of asymptotically free integrable 
theories should allow to understand how do the usual tools of modern theoretical 
physics (functional integration) apply to them. We consider this problem as the 
most important one: we should learn from exact solutions how to perform functional 
integration. There are also other reasons why such models as perturbations of 
WZNW model and principal chiral field are interesting for the application of form 
factor bootstrap which will be clear later. 

Suppose we deal with an massive integrable model which contains only one parti- 
cle in the spectrum (with isotopic degrees of freedom) with the two-particle S-matrix 
S(P) which satisfies Yang-Baxter, crossing, unitarity [4]. Then for the operator to 
be local it is necessary and sufficient that its form factors (matrix element between 
vacuum and n-particle state) satisfy the system of equations [1] which naturally 
splits into two parts: Riemann-Hilbert problem 

/(/?!,••• ,A,/W" ,/3„)S(A-A +1 ) = /(/?!,• •• ,&+!,&,-•• ,Pn), 

fifli,-" ,0n-i,0 n + 2m) = f(p n ,l3 1 ,-~ A-l), (1) 



i 



Typeset by AmS-T^L 



and residue condition 



2ni res / 3„ =(3 „_ 1+7 ri f(0i,--- , Pn-2, Pn-l, 

= /(/?!,••• ,P„-2)<S>S n - 1<n (l - S(p n -! 



^■■■SG^-l-A,^)) 



(2) 



where the usual [1] conventions are made: / belongs to the tensor product of spaces 
h* related to particles (h is one particle isotopic space), S(/% — /3j) acts nontrivially 
only in the tensor product of spaces related to particles with rapidities /3j , 
permutation of spaces is applied if corresponding rapidities are permuted, s n -i >n 
is a vector in the tensor product of (n — l)-th and n-th spaces constructed from 
charge conjugation matrix, / is meromorphic function of all its arguments in finite 
part of complex plane, as function of (3 n it does not have other singularities for 
< Im/3„ < 2iti but simple poles at the points (3 n = (3j + wi, j < n — 1. 

Several solutions to this infinite system of equations are known for the SG model 
which correspond to the most important local operators. We shall consider SG only 
in repulsive case when two component soliton is the only particle in the spectrum, 
so the space h is two-dimensional, and 4x4 S-matrix S^ G (/3) depends on the 
coupling constant £ : 7r < £ < oo. The explicit formula for this S-matrix will be 
given later. When £ is rational the restriction in the space of states is possible 
[5] which provides different model (4>i ^-perturbation of minimal model [6]). These 
restrictions do not respect the Hermitian structure of the space of states of SG 
and allow their own definite Hermitian structure for - integer. The S-matrix for 
the restricted models is obtained from SG one by RSOS [7] procedure. Generally, 
the idea of using RSOS restrictions for constructing physical S-matrices is due to 
[8] . In this paper we shall not go into much details of the restrictions because the 
restrictions of form factors occur quite naturally and do not present much difficulty 
as far as SG form factors are known. It should be mentioned also that for £ = oo 
the SG S-matrix produces the S-matrix of S'[/(2)-invariant Thirring model: 



There are several models for which the one-particle isotopic spaces are tensor 
product of two SG isotopic spaces (or their restrictions) the S-matrices are different 
particular cases (and restrictions) of the following one 



for two different coupling constants [9,10,11,12]. The most interesting examples are 
the following. For £i = 00,^2 = k + 2 after restriction of the second S-matrix we 
deal with the perturbation of WZNW-model on level k with action 



which is the same as fc-flavour 5?7(2)-Thirring model due to Polyakov-Wiegmann 
bosonisation [9]. In extreme case £1 = 00, £2 =00 we get principal chiral field model 
(PCF) [9,12]. These are the models we are mostly interested in. To get the form 
factors for them we shall consider the most general S-matrix of the type given by 




(/?) 



(3) 




(3). 
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Trying to solve the equations (1,2) for the S-matrix (3) we find the following 
amusing circumstance. Consider the Ricmann-Hilbcrt problem (1) for the S-matrix 
(3). Evidently, it is satisfied by a slight modification of the tensor product of SG 
form factors: 

exp(l^ft) Hcth^i - (3 j)fM,--- ,ft) (4) 

i<j 

Also this function is meromorphic and does not have other singularities as function 
of ft but usual simple poles at ft — (3j+iri (SG form factors vanish when ft = f3j). 
However, this anzatz breaks the equation (2). So, we have to look for something 
more intelligent. The lesson we learn from the naive anzatz (4) is that if we consider 
any solutions to the Riemann-Hilbert problem (1) for SG S-matrices then their 
tensor product satisfies the same equations for the tensor product S-matrix, but we 
should take care of the third equation. 

On the other hand the Riemann-Hilbert problem (1) for SG model can be consid- 
ered as deformation of Knizhnik-Zamolodchikov [13] (KZ) equations on level zero 
for the algebra U q {sl{2)) with q = exp(^- i ). This is rather informal way of think- 
ing because strictly speaking the deformations of KZ equations (thorough vertex 
operators, highest weight representations etc) are properly defined for q < 1 [14], 
for Yangian case q — ► 1 they can be treated in terms of asymptotic series [15], for 
\q\ = 1 we still do not know how to treat them, but fortunately on level zero we do 
know how to solve them! 

The same equations for the S-matrix (3) can be called deformed KZ for non 
semi-simple algebra U qi (sl(2)) ® U q2 (sl(2j) with q n = cxp(^- L ). Certainly the 
solutions for this case should be given by linear combinations of tensor products of 
different solutions of equations for U qi (sl(2)) and U q2 (sl(2)). So, our goal will be 
achieved if we know enough solutions in these two cases in order to construct linear 
combination of their tensor products which satisfy the residue condition (2). As it 
had been mentioned in [1] and explained in details for Yangian case in [15,16] we do 
know many solutions to the Riemann-Hilbert problem (1) in SG case (actually the 
same number as for usual KZ equations), but only very special ones were used for 
SG form factors because we had to satisfy the residue equation. In this paper we 
shall show that all the solutions are needed in order to construct the form factors 
for the S-matrix (3) through the procedure explained above. 

2. Solution of Riemann-Hilbert problem for U q (sl(2)), \q\ = 1. 

In this section we shall be interested in the solutions to the Riemann-Hilbert 
problem (1). It is convenient to change the sign in the RHS of the second equation, 
this is harmless because the solutions to modified in this way equations can be 
transformed into the solution of original ones via multiplication by exp(i^ft). 
So, we want to solve the equations: 

/(ft,--- ,&,&+!,■■■ ,02 n )S 6 (0i-0i+l) = f(pl,--- ,ft+l,ft,--- ,ft„), 

/(ft,--- ,ft„-l,ft„ +2ni) = -/(ft„,ft,--- ,ft„-l) (5) 
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where the SG S-matrix is given by [4] 



oo . 

qm , %o(/3) g , a , v (ft , , . /- sin(fcg)sh(^fc) 



^,o(/3) 
sh|(/?-7n) 



/ sh|(/3-7ri), \ 

-sh|/3 sh^ 

sh^ -sh|/3 

\ shf(/3-7ri)7 



Later we shall omit index £ when only one U q (sl(2)) is involved. We shall consider 
the solution to these equations of special isotopic character. The algebra U q (sl(2)) 
contains two finite dimensional subalgebras isomorphic to U q (sl(2)). The isotopic 
spaces of particles can be considered as the spaces of two-dimensional representa- 
tions for these two subalgebras (the explicit formulae can be found in [5] ) . We shall 
restrict ourselves with the consideration of those solutions which belong to invariant 
with respect to one of these subgroups spaces in the tensor product (the solutions 
will carry an index 7 = ± to indicate that they are invariant with respect to one or 
another subalgebra). For SG case these solutions describe the most fundamental 
local operators: energy- momentum tensor and current. This consideration can be 
generalized for \q\ < 1 case [17]. 

Let us mention that the solutions to the equations (5) are defined up to multi- 
plication by arbitrary quasiconstant: 27ri-periodic, symmetric function of [3j. The 
solutions which will be presented later are supposed to constitute the full set of 
meromorphic solutions of given isotopic structure up to quasiconstants. 

For what follows we shall need two special functions: <p(/3) and C(P)- We shall 
not write down explicit formulae for these functions which can be found in [1], but 
just present their most important properties: 

sh-(3 - 1 
- 2tt») = - * * v (/3-mM0) = 



shf(/3-^f)' rv ' rv ' 2sh|(/3-f)sh(/3-f ) 
if^t] = W), C(P)SoW = C(-/?), 

HTO -I 

C(/?-27ri) = C(-/?), C(/?)C(0-7ri) = (¥>(/? +y)) (6) 

We shall also need the vector-functions ^( a b ' ' ' , • • • ,p2n) defined by 

the following list of requirements: 

1) they take values in the tensor product of 2n spaces h* (h ~ C 2 ), 

2) they are antisymmetric, entire, periodic with period 2£i functions of cti , ■ ■ ■ , a„_i 

3) they are antiperiodic with period 2£i functions of (3\ , ■ ■ ■ , (3 2n satisfying the fol- 
lowing symmetry property: 

32 • • • ,a„_i|/?i, • • • ,(3i,(3 i+1 , ■ ■ ■ ,/3 2 «)S'5(/3i - A+i) = 
= 32 ("l, ' ' • , a n-i\Pi, ■ ■ ■ ■ ■ ■ ,/hn), 

4) for 7 = + or — they are singlets with respect to one or another finite-dimensional 
quantum group as explained above, 
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5) they do not have other singularities but simple poles at the points (3j = f3i + iri + 
k^i for j > i, fceZ 

6) they satisfy the recurrent relations: 

reS/3 2?l=/ 3 2rl _ 1+7r i dZ( a ir-- > <Xn-l\01, • ■ ■ , p2n-2, p2n-l, fan) = 

= —(e2n-l - ® e 2n .+ + e 2 „-l,+ ® €2n,-)® 
7T 

n-1 

(X]( _1 ) i ^-i( Q!l ' ,a«-i|/3i, • • • ,/3 2 „- 2 ) exp(-(n - l)(aj - ftn-i)) 

;=i ^ 



2n-2 . 2n-2 

X 

P 

n-1 . 2/i 2 

X 



[ n ^ n sh >< - & - ?)] ) 

=i ^ ^ P =i ? ^ 



FT Shj(a g - /3 2n _x - ?)( TT ^T^n-l - Pi)) 1 (7) 
g=l * ^ i=l * 

where e^± is basis in i-th space, the vector 

e 2 „-i,- <8> e 2 „,+ + e 2 „-i,+ (g) e 2 „,- (8) 

is nothing but s 2n _i 2n from (2). 
7) for n=l we have 

37(/Ji, ft) = 

cxp(7^(/3 2 - /?i - 7ri))(ei i+ ® e 2 _) + exp(~7^(/3 2 - fa - 7ri))(ei _ <g> e 2 , + ) 

sh£(/? 2 -Tri) 

It is possible to give explicit formulae for these functions in terms of certain 
determinants [1], we shall need it only in the last two sections of this paper. 

Now we are ready to write down the solutions. Different solution will be counted 
by 7 = ± which specifies the isotopic structure (singlet with respect to one or 
another finite-dimensional quantum subalgebra) and a set of integers {k\, • • • , fc„_i} 
such that \ki\ < n — 1, Mi. The sets {k\, • • • , k n -i} play the same role as different 
contours in integral formulae for solutions of KZ equations [13,18,19]. The solutions 
are given by 

fi u ...,k n - 1 (h,~-,02n) = <r UWi-b) (9) 

i<j 

+ ?° + ?° n-1 2n n-1 

/ dax--- / da„_i JJ JJ<p(ai-/3j)exp(^aiA;i)5^(ai,--- ,a„_i|^i,--- ,/3 2 n) 

J J i — 1 A — 1 A — 1 



i=l j=l 



where ^ 

The integrals in (9) are regularized in a special way explained in [1], this regu- 
larization makes sense only for hi integer, the limitation \ki\ < n — 1 comes from 
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requirement of convergency of regularized integrals. As it is shown in [1] the func- 
tions ;fcri _ 1 • • • ,/?2n) do satisfy (5). 

Let us make some comments on the sets {fci, • • • , k n -i} which count the solu- 
tions. First, due to antisymmetry of §1 with respect to a's we must consider only 
those sets with all fc, different. Second, for "good" entire, 2£i-periodic function 
F(a) we have the following relation: 



n +X° 2n 

X^i-i^V" .e*") / da JJ <p(a - /?j)exp((n + 1 - 2i)a)F{c 

»=1 „ .7=1 



:) = (10) 



where <Jk is elementary symmetrical polynomial of degree k. It is explained in [1,16] 
what kind of functions is "good", we would not go into details here, at least the 
function considered as a function of one a is "good". The formula (10) gives 
one (and only one) relation of linear dependence with quasiconstant coefficients 
between the integrals with different k. It provides, certainly, linear dependence 
between some solutions to (5), to get really different solutions we could, for example, 
require fcj ^= n — 1 \fi. 

There is a beautiful way of understanding the solutions (9) explained in details 
(for Yangian case) in [15,16]. The point is that the solutions of KZ equations for the 
case of sl(2), on level zero are given by determinants composed of periods of certain 
second kind differentials on hypcrclliptic surface. The size of these determinants is 
equal to genus of the surface while the number of different cycles in twice bigger, 
that gives rise to different solutions counted by different subsets of cycles. The 
deformation of this picture should be understood as follows. The points /3i, • • • $2n 
are branching points of "quantum hyperelliptic surface", the integral 

+ °° 2n 

ex P(-jE&)*2i-i(e ft ,---,e A ") / da f[ ip{a - (3j)exp{ka)F(a) 

is the period of differential defined by 2£i-periodic function F(a) taken over the 
cycle around two branching points: (3 n -k, Pn-k+i- The "good" functions men- 
tioned above are analogues of those differentials which do not have simple poles 
on surface. If we consider first the Yangian limit (£ — > oo) and then classical 
limit of Yangian then the correspondence can be explicitly understood in terms of 
asymptotics [15,16,20]. 

3. Form factors for SG model and its restrictions. 

Let us return to SG model. The form factors should satisfy not only the 
Ricmann-Hilbert problem (1), but also the residue condition 

2ni reS / 3 2 „ =/ 3 2 „_ 1+Tr i • • • , 02n-2, 02n-l, (hn) = 

= f(01, ■ ■ ■ A-2) ® *2n-l,2„(J - S(0 2n -1 ~ Pi) ' ■ ■ S(0 2n -1 ~ /?2n- 2 )) 

So, we have to consider the problem of calculation of residues of the kind. 

It is easy to show that generally the expressions (9) have simple poles at fan = 
/fen-i + 7ri. The recurrent relations (7) are useful for the calculation of residues 
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because together with the properties of ip (6) they provide that the integral with 
respect to ai in l-th term of (7) substituted to (9) can be replaced by contour integral 
over the contour: (-co, oo) (oo, oo + 2iri) (oo + 2ni, — oo + 2ni) (— oo + 2ni, -co). 
So, this integral is calculated via the poles of the integrand which are situated at 
the points ai — fii n -\ + a i = ftn-i + ^f 1 - However, generally the expressions 
for the residues can not be expressed in terms of functions of the type (9) with 
2n — 2 points f3j. There are only three possibilities to combine the functions (9) 
with different {fa} in order to have nice expressions: 

+ ?° + ?° n-l 2n 

2m res 02n=/32n _ 1+wi d n JJ C(ft - ft) / da x ■ ■ ■ j da n _ 1 JJ J| ^(a, - /?,•) 

'<i -oo -oo <=1J=1 

n— 1 n—1 

x J| 2ch(a 4 - /3 2 „-i)exp(^ a 4 ^)^(ai, • • • ,a„_i|/?i, • • • ,/3 2 „) = 0, 



+ 00 +00 



da!--- / rfa n _i J| J| yj(ai 

. ' I .5 = 1 



ft) 



n-l 



x (1 + ie" 1 "^" 1 ) J] 2ch(a 4 - /3 2 „_ 1 )exp(^ a % k) (11) 

i=2 i=l 

x S^(ai, • • • ,a„-i|/3i, • • • ,ft„) = 2s 2n _i ;2 „ ® fZ,-,i n S^ ' ' ' 'ft™)' 



+ ?° + ?° n-l 2n 

2tt« res / 3 2 „ =/ 3 2 „_ 1 +„ d" J|C(ft-ft) / d«i • ■ ■ / da n _i - ft) 

*<J -oo -oo 1=1 J = l 

n—1 n—1 

x (1 - ie" 1 "^"- 1 ) [| 2ch(a 4 - /3 2n _i)exp(^ a t k) 

i=2 i=l 

x S^(ai, • • • ,a„_i|/3i, • • • ,(3 2n ) = 

= -2S2 n -l, 2n ®f? li ... iln _ 1 ((h,--- ,ftn)S(ft«-l ~ Pi) • • • S{fo n -1 -ftn- 2 ), 

The limitations for the values of Zj are clear in every case. 

The solutions for the complete system of relations (1,2) for n > 1 are given by 



/+(ft, ■ • • , fan) = CXp(i ft)/-(n-l),-(n-3),... ,(n-5),(n-3) ' ' ' ' ' ft™)' 
/+(ft, ■ • • , ftn) = CX P (i ^ ft)/Z (n _3),-(„- 5 ),... )(n _ 3 ),(n-l)(ft ' • ' ' ' ft™) 



(12) 



To prove that these functions do satisfy the residue condition we notice that due 
to antisymmetry of ^ under the integral in (9) the expression 



n-l 



exp(^(n + 1 - 2i)ai) 



i=i 
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can be replaced by 

icxp((n - 2)ft n _i - (n - 3)ai) [(l + ie 01 ^ 2 "- 1 ) + (l - ie 011 -^"- 1 )] 

n— 1 n— 1 

x 2ch(ai - /3 2 „-i)exp(^(n - 2i)aj) 

i=2 i=2 

Now one just uses (11) to prove that the residue condition is satisfied. 

The functions (12) are not independent, due to (10) they satisfy the relation 

Analyzing the two particle form factor one realizes that the form factors of 
energy-momentum tensor T ai(72 and U(l) current J a (we use light-cone components, 
a = ±) are given by [1]: 

2n 

/++(/?!,••• ,/3 2 „) = eft) (ft,--. ,&„), 

J=l 7 
2n 

J=l 7 

/+-(&,••• ,A„) = 

2n 2n 

= (E e ~ ft ) E /+(a, ■ ■ ■ > a.) - (E eft ) E /-(a. ■ ■ ■ . An). 

j=l 7 j=l 7 

/±(/3i,--- ,A„) = E(-) 7 /2(/?i,-- - ,A„) (13) 

7 

There is one more set of solutions which have good residues: 

f-(n-2),-(n-4),- ,(n-4),(n-2)(A) 'An) ( 14 ) 

but in that case we keep minus in the RHS of (1), also we have plus instead of 
minus between two terms in the RHS of (2), hence these solutions corresponds to 
disorder operators [1]. 

Let us explain briefly the restrictions of SG. Consider the modified energy- 
momentum tensor whose form factors are given by 

2n 
2n 

/|-(A,"- ,An) = 

2n 2n 

= (E e_ft )/+(A,--- ,An) = (E e/33 ')/-(A.--- .An) 

j=i i=i 



These operators is invariant with respect to one of U q (sl(2)) subalgcbras. That 
is why for rational £ the intermediate states in the correlations of these operators 
among themselves happen to be truncated [5] . So, we can construct a new restricted 
theory (RSG) with smaller operator content and truncated space of states. This 
theory is well known to coincide with <j>\ ^-perturbation of minimal model [6]. The 
truncation does not respect the Hermitian structure of SG space of states, so, to 
have RSG to be defined intrinsically we have to introduce new Hermitian structure. 
It can be made positively defined in the region on coupling constant we deal with 
only for |- integer and £ = | (in what follows we shall not be interested in the latter 
case). One of disorder operators allows restriction, it coincides with the operator 
<f>i_2 in ultraviolet limit. 

Thus, only very small part of solutions we know is useful for SG form factors. 
However, in the next section we shall show that all of them are needed for the 
models with the tensor product S-matrices. 

4. The form factors for tensor product S-matrix. 
The Riemann-Hilbert problem (1) looks in that case as 

f(Pl, • ' • i Pi, Pi+1, • • ' ) 02n)S£i&(0i ~ Pi+l) = f{Pl, ' ' ' , Pi+1, Pi, ' ' ' , Pin), 

with 

S 5l5 M = Sl G (P)®Sl G (P) (15) 
Evidently, these two equations are satisfied by any expression of the kind: 

exp(±i &) II ctl 4(& - Mf£, *(&>•■■. &™) ® Ul l(&> ■ ■ ■ . &») 

i<j 

where we denoted the ordered subsets of integers {fci, • • • , : < n — 1 and 

{h, ■ ■ ■ ,l n -i} ■ \h\ < n ~ 1 by K and L. The problem is to satisfy the residue 
condition: 

2ni reS/3 2?l=; 3 2rl _ 1+7r i f{Pl, • ■ • , P2n-2, p2n-l, P2n) = 

= f(Pl, • • • )#2n-2) ® S2n-l,2n{l ~ <%?2 {02n-l ~ Pi) ' ' ' >%?2 (An-1 — /?2n-2)) 

(16) 

where S2 n -i,2n is the tensor product of two vectors like (8). We shall show that 
these equations are satisfied by 

/±66 (ft>"-> Pm) = (2^) n cxp(±i &) II ct 4 (ft ~ & 

i<j 

X E f£,K(Pl,--- ,p2n)®fg, L (Pl,--- ,p2n) (17) 
N±=KuL 

where N+ = {-(n - 1), -(n - 2), • ■ ■ , (n - 3), (n - 2)} , N_ = {-(n - 2), -(« - 
3), • • • , (n — 2), (n — 1)}. To prove that we have to calculate the residue. It is not 
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complicated, first let us write more explicit formulae: 

- ((n \ w d l^M±lT,h)Il^M ~ &) x 

+ ?° +?° +?° + ?° n-l 2n n-l 2n 

/ dai • • • / da n -i / dan ■ ■ ■ / dd„_i J j ]^[ </> €l (a* - ]J j J (d; - fy) 

-oo -oo -oo -oo i=1 J =1 i=1 J =1 

x exp(T^J2( ai + &i ^n 2sh \( ai ~ a i)II 2sh ^ 5!i ~ "j')II 2cl 4( ai ~ 

xt fl K- ,a„_i|/3i,--- ,^ n )®3^£ 2 ( 5 i.-" . 5 n -i|/9i," • ,/3 2 n) (18) 
where 

dfi& = ™* €l d &) = cthi(/3)C 6 (/3)fe(/3) 

Let us prove that for n > 1 the functions (17) do satisfy the residue condition (the 
case n = 1 is special). To do that it is sufficient to realize that under the integral 
due to antisymmetry of $ we can perform a replacement: 

n shi( a , - aj -) nsh^ - a,-) n n <4(«< & ^ - 2- 2 ("- 2 > e <»- 2 >*»-i 

i<j i<j l<i l<j 

1 1 

x {(-l)"(n - l) 2 cxp-(-(n - 2)(ai + d 1 ))ch-(a 1 - di) J] ch(a, - /3 2n _i) 

i=2 

n-l 111 

x J| ch(di - /fen-i) I| sh-(ai-aj) J| sh-(di - dj) ch-(a!; - dj)+ 

i=2 2<i<j 2<«<j 2<i2<j 

(19) 

(n — l)(n — 2) 1 , > 1 ^-r 1 

H ^ cx P2( _ ( n - 2 )( ai + a 2)jsh-(ai - a 2 ) [[ ch(a; - /3 2n -i) 

i=3 

n-l 111 

]J ch(di - fo n -i) Y[ sh-(ai-aj) ]| sh-(d 4 - dj) J| ch-(a, - dj)+ 



n-l 

X 



2 xx 2 xi ii 2 

j=l 3<i<j l<«<j 3<i l<j 



(n — l)(n — 2) 1 , _ > 1 Vr 1 

H ^ exp-(-(n - 2)(ax + a 2 ))sh-(ai - a 2 ) || ch(a; - /? 2 n-i) 

i=3 

n-l 111 

x JJch(aj-^2 n _i) Y[ sh-(d i -d J ) J| sh-( ai - aj ) ch-(d; - aj )} 

i=l 3<i<j l<i<j 3<il<j 

This formula looks quite terrible, but it has a simple origin: the LHS is proportional 
to Vandermonde determinant composed of e ai and — e ai , so, to get (19) we added 
in this determinant to fc-th (for k > 3) row the (fc — 2)-th one multiplied by e 2/32n_ 1 , 
and then decomposed with respect to first two rows. Notice now that we are 
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interested in the second order pole at ft„ = fiin-\ + of the integral from (18) 
since Cfi.&Oftn — ftn-i) has zero at this point. But the contributions from the last 
two terms in (19) do not have such singularity. Consider, for example the last term. 
The integral with respect to d's will produce first order pole, but the integral with 

n-1 

respect to a's is regular because it contains Yl ch(«i — /?2n-i) ( sec (H))- Thus 

»=i 

the only interesting term is the first one from (19). In this term we can replace 
chi(«i - di) by 

(1 - ze Ql - /32 - 1 )(l + ie- 31 ^ 2 "- 1 )] 

and then use the formulae (11). That proves the formula (16) for n > 2. Let 
us mention one more important property of these form factors: they satisfy the 
relations: 



2n 2n 

\ f7i72 I a ... a \ = AT P ftW7i72_ 



i=i i=i 



these relations are proven using the formula (10). As usual, to understand what 
kind of operators do these form factors correspond to we have to investigate better 
the two-particle ones. But before doing that we have to explain what kind of theory 
we deal with. The point is that the S-matrix in question describes many different 
models. 



5. Different models described by tensor product S-matrix 

The simplest way is to understand the S-matrix as it is, and to introduce the 
scalar product in the space of states such that it respects the Hermitian conjugation 
of S-matrix: 

Then the symmetry of the model is U(l) <g> U(l), two-particle form factors are given 

by 



= d 5l?2 exp(±i^/3,)-— 



-ft) 



sh£(ft - ft - ™)sh£(ft _ ft _ m ) 
® [ e -»*^- ft -^ (e li+ ® e 2 ,_) + e -72 ^(fe-ft— ) ( ei _ e2>+ ^ 



From these expressions an from the formula (20) one realizes that the form factors 
of energy- momentum tensor T ai(T2 and two U(l) currents J^, (we use light-cone 
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components, a — ±) are given by 



2n 

• • • , fo n ) - (E eft ) E /+ 66 0^1. • • • . 

j = l 7172 
2n 

/__ (ft, ■ ■ ■ , /3 2 „) = (E e- ft ) E /-66 (A . ' ' ' > M, 

j=i 7l72 

/+-G8i,"- ,&„) = 

2n 2n 

= (E e " ft ) E /+ 66 (A . • • • . A») = (E eft ) E /+ 66 (A. ■ ■ ■ . 

J=l 7172 i=l 7172 



J±(/?l, • • • ,02n) = E(-) 7l /±66^' • • • " 32 ")' 
7172 

/£(/?!,■■■ = E(-) 72 /r^ 2 (/3i,--- ,/3 2 „) 



7172 



The important limit of this model is £i,£2 ~~ * 00 ■ ^ n thi s limit we get the PCF. 

Now we can consider two levels of restriction of the model. First, let us introduce 
the modified energy-momentum tensor and modified left current as 

2n 

ft+ih, ■ ■ ■ , fan) = (E eft ) E /+ 66 ■ ■ ■ ' 0*O> 
j=l 7i 

2n 

/±_ (a, • • • , fa) = (E e ~ 0j ) E /-66 (0i. • • • < &») > 

j=l 71 

/|-(01,"- ,02n) = 

2n 2n 

- (E e " ft )E/+66(0i.--- .a») = (E eft )E/+66(0i.--- >&»)> 

j=l 71 j=l 71 

/±' + (/3l,-" ,02n) =EH 7l /^e 2 (01.-" -02"), 



These operators are not selfajoint in SG Hermitian structure, but they are invariant 
under the action of one U q2 (sl(2)). That is why for rational — the correlations of 
these operators among themselves contain only RSOS restricted with respect to 
Uq 2 {sl{2)) states. For | = H2,fceZ the model with restricted set of operators 
and truncated space can be equipped with positively defined scalar product. It 
is no wonder that is lost in the restricted model: right U(l) is broken. The 
important limit of the restricted model is £i — > oo for given k. In this limit we get 
perturbations of WZNW on level k. 

It is instructive to recover SG model itself. It should coincide with the restricted 
model for £i = £, £ 2 = 37T. In that case in the truncated space of states right 
degrees of freedom are frozen and the restriction of S^ 2 (P) is just —1 (Ising S- 
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matrix). Consider, for example, the formula 

/!+(/?!, • • • , An) = (2^)"(^eft)exp(±- II cth 2^ ~ M 

*E E /?, K (PU- L (A,"-,An) 
7 N + =K\JL 

It can be shown that the restriction of 

/^(A,-" ,An) 

differs from zero only for L = {n — 2j, j = 1, • • • , n — 1}. But in the latter case 
the restriction coincides (14) with the form factor of Ising disorder operator which 
is given by 

this expression cancels in (17) and we recover SG energy-momentum tensor form 
factors (13). 

At last if we consider another modification of the energy-momentum tensor with 
form factors 

2n 

/++(A, ■ ■ ■ , An) = (E eft )/tU (A, • • • ' An), 
2n 

• • • ,An) = (£ e ~ ft ')/- • ' ' .^«). 
/+±(/3 1 ,-- - ,/3 2 „) = 

2n 2n 

then for ^ , ^ integers we can perform restriction which will give the perturbations 
of coset models discussed in [11]. 

It should be said that we restricted ourselves in this paper with consideration 
of repulsive SG coupling constants (£ > it). The consideration of attractive case is 
also important, after restriction we can find quite unexpected models. 

6. The mathematical structure of the solution. 

The S-matrix —S^({3) <g> S'£ G (/3) is constructed in amusing way: as if we have 
two types of particles confined together. The formula for form factors 

/± (A, • • • , An) = exp(±i ]T A) II cth ^(A - A) 

i<3 

X E /€ 1 U(A,---,An)®4 2 , i (A,---,An) (21) 
N±=KuL 
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shows that these two types of particles "interact" in quite interesting fashion. In 
this section we shall explain the most attractive features of this " interaction" . 

When considering the solutions to Riemann-Hilbert problem we explained that 
the sets K which count different solutions play the same role as choice of different 
contours in integral formulae for the solutions of KZ equations. From the algebraic 
point of view these data count different passes composed of intermediate Verma 
modules for the products of vertex operators whose vacuum expectations provide 
the solutions to KZ and deformed KZ equations. This interaction through passes 
looks quite interesting, it is similar to combining together holomorphic and anti- 
holomorphic conformal blocks in CFT. 

However, we suppose the following circumstance to be of the main importance. 
We deal with k = equations, in that case the contours in question are quite 
special, as it has been mentioned in Section 2. Let us explain this point in more 
details, but before doing that we have to provide some more information. 

Let us concentrate on one SG S-matrix with coupling constant £. First, we shall 
give more explicit description of the vector 3^(ai, • • • , a„_i|/3i, • • • , /?2n)- 

We introduce the notation B = {1, 2, • • • , 2n}. To every T C B : j^T = n vector 
wt{(3\, ■ ■ ■ , 02n-i, @2n) from the tensor product of isotopic spaces is related such 
that 

W T (01, ■ ■ ■ ■ ■ ■ ,/?2n - = W T \i{01, ■ ■ ■ ,0i+l,(3i, ■ ■ ■ ,02n), 

where 

f T if i,i + 1 e T or i,i + 1 £ T 

Toio{i + l} if i e T, i + 1 £ T or i + 1 e T, i <£ T 

where AoB = (AL)B)\(Ar)B). These notation are different from those used in [1], 
but one can easily understand the relation. It would simplify a lot further formulae 
if we introduce notations 

We have [1]: 

$Z(ai, ■ ■ ■ ,a„_i|/3i, • • • ,/3 2n ) = 

= A T («!,-•• ,a„_i|/3i, ••• ,/3 2n )fJ(/3i, ••• ,/?2n) 

T£B,#T=n 

where A T is (n — 1) x (n — 1) matrix with matrix elements: 

Aj, = AfiajlPu--- ,/3 2n ), 



where 



Af (a|/3 1; • • • , (3 2n ) = exp(-|((n - 2)a + ^/3 P )) 



i-l 



q£T k=0 



+ q i H ( e T« - g-lerA) ^(1 - ^JeT^^ 1 ) ,,!^ 

qeB\T k=0 
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with <7fc ; T, Cfe, b\t are elementary symmetric polynomials of degree k with argu- 
ments exp^ L (^/3 p ), p eT and exp^ L (J]/3p), p respectively. We remind also 
that q = exp^i. 

These formulae provide that the functions k(0i> ' ' ' iftn) can be written in 
the form: 



wj(/3i, • • • ,p 2n ) 
(22) 

where K — {hi, • • • , fc„_i}, fci < /c 2 < • • • < k n -\. 

To rewrite this answer in more beautiful way let us introduce two vector spaces: 
the space H of dimension In — 2 with basis Ai, Bi i = 1, • • • , n — 1, and the space V 
of dimension n — 1 with basis Zi, i = 1, • • ■ , n — 1. Now for every T C B, #T = n 
we introduce two forms 

w ± (T) = u id (T)Ai A Zj ± Uij (T)Bi A Z,- (23) 

where 

n— 1 n— 1 

WiJ (T) = cxp(-- ^ /3 P ){£ ^ <7 2 „_ 2j (e^, • • • , e*-) 

i=i j=i 

x / JJ<^(a-/3 p )Af(a|/3i,--- , /? 2 „)exp((2p - n)a)da} 

(T) = exp (-i &){£ E E ct2 - 2 '+! . • ■ ■ ' e/3a ") 

i=l j=l 1=1 

x / n^-^^HAr-- ,/3 2 „)exp((2/-n-l)a)rfa} 
J p 

Then different solutions (22) can be found as coefficients in decomposition with 
respect to 

(±) q A h A • • • A A jq A B kl A • • • A B kn _ q _, (24) 

of the expression 

]T A^^mKGSi,--. ,/3 2 „) (25) 

T£B,#T=n 

We can think of this decomposition as of one with respect to Grassmanian of (n — 1 )- 
dimensional subspaces of H. 

The "interaction" in the formula for the form factor /+ 1 ^^ 2 (21) can be under- 
stood as follows. Take two copies of the spase V (with basises Zi, Wi). Then the 
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E det 



jY[ip(a- (3 p )Af (a|/?i, • • • ,/3 2 „)exp(/cja)cfo! 



form factors are defined by the inner product of forms A^ n and A^ n 1 ^^ 2 as 

follows 

d iti ex p(f e ^) n (a - ^ ) ( n ^ • • ■ • - e/?2 - )) _2 x 

i<j p=l 

TieB,#T=nT 2 GS,#T=n 

A ( " _1) (W£ a , ij(T 2 )Ai A Wj- - %, ij(T 2 )Bi A Wj) 

X W^ iTi (ft , • • • , ft„) <g> «g iT2 (ft , • • • , (3 2n ) = 

xA 1 A • • • A A„_! A Si A • • • A B„_i A Zi A • • • A Z n _ x A W\ A • • • A W n _i 

Certainly, the formula (25) presents rather formal way of rewriting (22), but it is 
not an empty exercise since the form u)(T) independently has interesting meaning. 
It has been said in Section 2 that the expression 

exp(-^/3 P K_ fe (e^,-- - ,e*») 

x / JJ^(a-/3 p )Af(a|/3i,-- - , /3 2 „)exp(/ca)c?a 
^ p 

can be considered as quantum deformation of the period of special second type 
hyperelliptic differential over the contour around two Ck around two branching 
points (i n +k and /3 n+ k+i- The canonical choice of homology basis on the surface is 
the following: 

a-i = C2i- n , bi = C_„ + i + C- n+ 3 + • • • + C_„ +2 i-i 

for i = 1, • • • , (n — 1), genus of surface equals n — 1. From that point of view it 
is natural to identify Ai , Bi with basic vectors of the lattice of periods and the 
vectors Zi with differentials along Jacobian. Then the form u>(T) has nice meaning 
[21]. We shall explain this point in some more details in the next section. The 
combination of lo + with oj~ reminds one more time combining holomorphic and 
antiholomorphic pieces in CFT: the 6-cycles are imaginary ones. 

7. Remarks on the classical limit. 

Let us discuss the perturbation of WZNW model on level k. In that case we have 
to put £i = oo, £2 = k + 2, to fix 72 (say 72 = + ) and to perform restriction. Let 
us consider, for example, the form factors of one component of energy-momentum 
tensor: 

/ ++ (ft, • • • , ft„) = (2^)"cxp(l ]T ft) (]T eft) I] cthi(ft - ft) 

i<j 

x E Mft,-" ,ft„)®/!(ft,-'- ,ftn) (26) 
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where the notations are used 



' ' ' iftn) = (/ifc+2), L^ 1 ' ' ' ' ' ^2n)) restricted 



(27) 



Two pieces of (26) looking quite similar have very different analytical structure. 
The function /ic(ft,--- , ft n ) is a solution to Yangian Riemann-Hilbcrt problem 
i.e. that with rational in (3 S-matrix. That is why it makes perfect sense to consider 
the following classical limit for this piece [20] : 



P.i 



2ir\j 



In that limit the asymptotics of fx (Pi, ■ ■ • , Pin) is described by a solution to usual 
KZ equations on level zero. On the other hand the S-matrix for /|(ft, • • • ,p2n) 
depends typically on exp(-j^^), and the limit (27) does not make much sense for this 
function. Our dream would be to present the form factor as functional integral (with 
special boundary conditions) of the original action of the theory. The theory has 
two main features: it is asymptotically free and it contains WZNW term in action. 
These two features should lead to different effects: the asymptotic freedom should 
provide reasonable perturbation theory with respect to Plank constant while the 
WZNW term should provide nontrivial nonperturbative effects. The exact solution 
(26) shows that these two effects are combined in rather special way. Notice that 
the first piece of (26) is independent of k, also it allows the classical limit (27), H 
should be identified with Plank constant: it rescales the rapidities (logarithms of 
momenta) which makes perfect sense in asymptotically free situation. The second 
piece is the one depending upon k, it should be related to nonperturbative effects 
due to WZNW term. The formula (26) and the Grassmanian of the previous 
section suggest that by introducing certain fermionic field we should be able to 
treat these two pieces independently and perform the averaging over this fermionic 
field (summation over K, L) afterwards. We hope to return to the consideration 
of functional integral in feature, but now let us concentrate on the quasiclassical 
theory of f K (ft , • • • ,ft„). 

What we shall do now is certain extension of the consideration of papers [20,21]. 
Let us take the form u>(T) from the previous section for the Yangian case (£ = oo). 
Now we calculate the asymptotics (27) the result being [20] : 



n— 1 n— 1 



± (T) = 5252([t?A j AZ i ±[cT B jAZi 

i=l .7=1 / u 



(28) 



where eij, 6, are basic contours on the hyperelliptic surface t 2 = P(X) = JJ(X — X p ) 
of genus g = n — 1 , Qf is the following second kind differential: 



Si 



x { U( x - x p) 

i per 



+ l[(x-x P ) 

+ p$T 
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d U P e T ( x ~ X p) 
dX X n ^ i 



where [] + means that only polynomial part of the expression in brackets is taken. 
When substituting w ± (T) into (25) we can perform certain transformations. Little 
modification of formulae from [20] provides that (25) can be rewritten in the limit 
as 

StC' 3 A ( "- 1) (c2; ± (T)) ^t] (0) 4 £ T 

TeB,#T=n 

where Et is basic vector in the tensor product of isotopic spaces. We introduced 
Riemann thcta- function taking dj, bi for the canonical basis, rjx is even nonsingular 
(such that 6[i]t] (0) ^ 0) half-period related to the subset of branching points defined 
byT[22], 

C = l[(\ i -\ j f*A, (29) 

i<j 

A is given by the determinant 

A^ 1 



A = dot 



The form Cj (T) is rewritten as 



/ 



(n-1) X (n-1) 



£>(T) = J2 ^(^9 J log0M(O)A J Adz,±^(r! 4 ^ ; 5 i log0[7?T](O) + 2^ J )B i Adz J ) 

i=l j=l 1=1 

(30) 

where we replaced quite formally Zi by differential along Jacobian dzi. 

Let us consider the form on the Jacobian which interpolates between (30) for 
different T: 

n—1 n—1 n—1 

c£) ± (z) = ^2^2(d l d j \og9(z)A i A dzj ± ^2(Qi,ididjlog9(z) + 2m5 i j)B l A dzj) 

i=l j=l 1=1 

where z varies over Jacobian: z G C"" 1 /^" -1 x ^ n ^ 1 The mathematical meaning 
of this form can be explained as averaging of the simplcctic form induced on Jaco- 
bian after embedding into projective space by means of theta-functions of second 
order [21]. 

Probaaly the best possible situation takes place for PCF. Here both pieces of the 
form factor allow classical limit, and the form factors are special values (at even 
non-singular half-periods) of the form: 

h^C- & e(zfe{wf{^ n -^ lu+(z)) A (A^- 1 ' u>-(w)) (31) 

It is interesting that all the formulae of this paper can be considered as deformations 
of (31). 

8. Conclusion. 

To conclude this paper let us formulate several problems which, to our mind, are 
worth investigation. 
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1 . It would be interesting to generalize the considerations of this paper to other Lie 
algebras. The necessary preliminary information for SU(N) case can be found 
m[l] ; 

2. It is interesting to consider the deformation of the construction of this paper 
to the lattice models [23]. In that case different choice of contours in KZ is 
replaced not by introducing the exponents under the integrals, as it happened 
in the situation of this paper, but to introducing different theta-functions. The 
generalization for the lattice analog of PCF should be not complicated to find, 
however, the generalization for lattice version of perturbed WZNW on level k 
which coincides with intcgrablc version of 6-vertex model of spin 2k ,\q\ < 1 is 
not easy to describe because it will require the knowledge of solutions to elliptic 
version of Ricmann-Hilbert problem (1) which are not known (the same is needed 
for form factors of 8- vertex model). 

3. We suppose that the most important question is that of understanding the origin 
of form factor formulae in terms of functional integral. The results of this paper 
should be important in understanding of this problem. Formula (31) should be 
crucial for this goal. 

Acknowledgement I am grateful to N. Reshetikhin for interesting and stimu- 
lating discussions. 
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